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ABSTRACT 

• We report the detection of Voigt spectral line profiles of radio recombination lines (RRLs) 

pP toward Sagittarius B2(N) with the 100-m Green Bank Telescope (GBT). At radio wavelengths, 

astronomical spectra are highly populated with RRLs, which serve as ideal probes of the phys- 
ical conditions in molecular cloud complexes. An analysis of the Una lines presented herein 
shows that RRLs of higher principal quantum number (n > 90) are generally divergent from 
their expected Gaussian profiles and, moreover, are well described by their respective Voigt 
profiles. This is in agreement with the theory that spectral lines experience pressure broadening 
■ as a result of electron collisions at lower radio frequencies. Given the inherent technical dif- 

ficulties regarding the detection and profiling of true RRL wing spans and shapes, it is crucial 
that the observing instrumentation produce flat baselines as well as high sensitivity, high reso- 
lution data. The GBT has demonstrated its capabilities regarding all of these aspects, and we 
believe that future observations of RRL emission via the GBT will be crucial towards advancing 
our knowledge of the larger-scale extended structures of ionized gas in the interstellar medium 
(ISM). 
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1. Introduction 

Sagittarius B2(N) is a region undergoing active star formation that hosts a number of environments, 
including expanding H2O masers (Reid et al. 1988); regions of large complex molecules (Snyder et al. 
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1995); in-falling, outflowing, and rotating gas (Qin et al. 2008); and ultra-compact and extended HII regions 
(Chaisson 1973; Mehringer et al. 1993; Miao et al. 1995). Of particular interest are the radio recombination 
lines (RRLs) associated with these HII regions which, due to their optically thin radiation throughout the 
Galactic disk, may serve as accurate diagnostic probes of ionized gas at great distances. For example, if 
a source exists in local thermodynamic equilibrium (LTE) and the impacts from electrons are negligible, 
then the line-to-continuum intensity ratio of the RRL may be used to determine the electron temperature 
(Shaver 1980). Since conservation of momentum dictates that electron motions dominate the dynamics of 
an ionized gas, the electron temperature is considered to be a close approximation to the HII region's kinetic 
temperature, which may further be used to constrain the boundaries of other physical properties such as 
source size, flux density, rms electron density, and the magnitude of turbulence within the gas (Brown et al. 
1978; Roelfsema & Goss 1992; Rohlfs & Wilson 1996; Gordon & Sorochenko 2002). Through a gradient 
in electron temperature over Galactocentric radius, RRLs have been used to determine the trend in Galactic 
metallicity (e.g., Churchwell & Walmsley 1975; Wink et al. 1983; Shaver et al. 1983; Quireza et al. 2006). 

The shape of the spectral line can also provide important information about the dynamics of the ionizing 
gas which constitutes the HII region. The electrons of an HII region in LTE are expected to experience 
a Maxwell-Boltzmann velocity distribution and subsequently produce Gaussian line profiles via Doppler 
broadening. Indeed, most published RRL profiles are well-described by Gaussian profiles. However, for 
RRLs with large principal quantum numbers (n) arising from high density gas, pressure broadening from 
electron impacts is expected to be important (Griem 1967). The profile generated by pressure broadening has 
a Lorentzian curve, which has higher intensity in the line wings than a Gaussian profile does. In many cases, 
RRLs reflect both thermal (Doppler) and pressure broadening in the ionized gas. Consequently, the observed 
shape of such RRLs is a convolution of the individual Gaussian and Lorentzian components-a Voigt profile. 
Unfortunately, the detection of Voigt profiles from RRLs is often challenging. The Lorentzian component 
is expected to become important at higher n, and thus lower frequencies. Because the RRL intensity in this 
regime is comparatively weaker than it is at higher frequencies, the Lorentzian-broadened wings are often 
difficult to distinguish from instrumental baseline structure (e.g., Smirnov et al. 1984). Furthermore, it is 
also easy to distort the RRL line shape using data reduction techniques (e.g., Bell et al. 2000). 

While many authors have reported detection of Voigt profiles from RRL emission toward HII regions 
(e.g., Simpson 1973; Smirnov et al. 1984; Foster et al. 2007), few have obtained the high spectral resolution, 
high sensitivity observational data necessary to achieve unequivocal results. The difficulty in obtaining 
such data is due to several factors: (i) low signal-to-noise ratios (S/N) for which the enhanced line wings 
of the Voigt profile cannot be detected; (ii) poor spectral resolution, which prevents sufficient sampling 
of the RRL profile; (iii) the presence of nearby RRLs (e.g., Heo-, as seen in Figs. 1-4), which produce 
confusion within the spectrum; (iv) the broad beams of many telescopes (especially single-dish ones), which 
dilute and mix regions dominated by Doppler broadening with those which contain significant pressure 
broadened components; (v) large-scale motions, such as infalling and outflowing ionized gas, which produce 
asymmetric wings to appear on a thermal (Gaussian) profile; and (vi) poor spectral baselines. The last factor 
has arguably been the most difficult to circumvent as instrumental effects can distort the spectral baselines, 
making it very difficult to distinguish the line wings from the baseline (e.g., Bania et al. 1997). Nevertheless, 
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there is ample evidence that pressure broadening is important in HII regions. In particular, non-LTE radiative 
transfer models can only match the observed RRL intensity or line width for some HII regions if pressure 
broadening is included (e.g., Churchwell 1971; Balser et al. 1999; Keto et al. 2008). 

To our knowledge, the only unambiguous detection of RRL Voigt profiles are at low radio frequencies 
(<100 MHz) where the lines are detected in absorption toward a bright continuum object. Kantharia et 
al. (1998) detected Voigt profiles of carbon RRLs toward Cas A with an effective integration time of 400 
hours. More recently, Stepkin et al. (2007) detected carbon RRLs toward Cas A from the largest bound 
atoms in space with n ~ 1009 and an integration time of 504 hours. These studies were able to detect 
Voigt profiles for several reasons: (i) the very long integration times provided the necessary sensitivity; (ii) 
the Lorentzian component of the Voigt profile is more pronounced at the lower frequencies observed; and 
(iii) the bright background intensity from Cas A and the Galactic background produces radiation broadening, 
which also contributes to the Lorentzian component. Measuring the Voigt profile was critical in constraining 
the physical conditions in these regions, which are thought to be associated with the cold atomic hydrogen 
component of the diffuse ISM. 

The design of the GBT, which leaves its aperture unblocked, provides many improvements in sensitiv- 
ity, spectral resolution, and spectral baselines over traditional single-dish telescopes. For example, multi- 
path reflections from the superstructure are avoided by the GBT (Fisher et al. 2003). Although much of the 
instrumental baseline structure produced by radio telescopes is correlated out in radio interferometers, until 
recently most instruments did not have sufficient spectral resolution to resolve Voigt profiles. 

In this article we report on the clear detection of (at least) four emission line RRL Voigt profiles from 
observations directed toward Sagittarius B2(N) as part of the Prebiotic Interstellar Molecule Survey (PRI- 
MOsQ). A sample of the 12 highest S/N RRLs detected in this survey are presented, with frequencies ranging 
from 4.7 to 24.5 GHz, corresponding to principal quantum number n between 64 and 111. All told, there 
will be about 1200 H RRLs present within the PRIMOS survey frequency range of 1-50 GHz. This includes 
only those RRLs from the An=l (alpha lines) to An=6 (zeta lines) series. The RRL rest frequencies were 
taken from the SplatalogueH database (Remijan et al. 2007). It is beyond the scope of this paper to discuss 
all the currently detected RRLs toward the Sagittarius B2(N) region. Rather, the main point of this work 
is to present the data repository where researchers can obtain the full complement of RRLs detected and to 
illustrate that Voigt spectral-line profiles of RRLs can clearly be detected with the GBT. As such, we discuss 
only a subset of the data available and illustrate that high-rc (n >90) H RRLs toward Sagittarius B2(N) are 
best fit with Voigt spectral line profiles. 

Figures 1-4 and Tables 1-2 summarize our results. Spectra of 12 Una RRLs are shown with Gaussian, 
Lorentzian, and Voigt function profile fits to the raw data. An expanded view of the line wings near the 
baseline is shown to illustrate the differences in these models. All spectra are plotted with respect to the 
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local standard of rest (LSR) velocity, centered at km s . Table 1 lists the fitting parameters for each of 
the Una RRLs. We present a review of the observations made in §2, describe our reduction methods and 
analysis in §3, and discuss the data and profile fits in §4. 

2. Observations 

All observations were obtained between 4.7 and 24.5 GHz with the NRAC^I 100-m Robert C. Byrd 
Green Bank Telescope from 2004 March 4 to 2005 November 12. Briefly, the data were obtained using the 
GBT spectrometer, configured to provide four intermediate frequency spectral windows simultaneously in 
two orthogonal polarizations through the use of offset oscillators. The spectral resolution was 24.4 kHz at 
frequencies centered between 4.7-24.5 GHz. This is significantly better than many previous measurement^ 

The GBT's half-power beam-width can be approximated by: 9b ~ where the observing frequency 
(y) is measured in GHz. Over the course of the observations, the spectrometer bandpasses were centered at 
different frequencies in the range of 4.7 GHz to 24.5 GHz that correspond to 9b ~ 157" and 9b ~ 30.2", 
and a beam efficiency from 0.96 to 0.83, respectively. The aperture efficiency is 0.70 to 0.61 over the same 
frequency rang^E 

The J2000 pointing position employed was a = 17 /, 47 m 19 s . 8, 5 = -28°22' 17" {€ = 0.677° , b = -0.027°), 
which is centered on the Sagittarius B2(N-LMH) source and an LSR source velocity of +64 km s _1 was 
assumed. The dynamic pointing and focusing corrections were automatically updated and employed based 
on realtime temperature measurements of the structure input to a thermal model of the GBT; local pointing 
and focus corrections were adjusted typically every two hours using standard pointing sources. Spectrometer 
data were taken in the OFF-ON position-switching mode, with the OFF position at 60' East in azimuth with 
respect to the ON source position. A single scan consisted of two minutes in the OFF source position 
followed by two minutes in the ON source position. Sagittarius B2(N) was observed in this manner above 
10° elevation from source rise to source set (a six-hour track). When available, all the accumulated scans 
of the two polarization outputs from the spectrometer were averaged in the final data-reduction process to 
improve the S/N. The antenna temperatures of the RRLs presented in this paper are on the T^ scale (Ulich 
& Haas 1976) with estimated 20% uncertainties. 



3 The National Radio Astronomy Observatory is a facility of the National Science Foundation, operated under cooperative 
agreement by Associated Universities, Inc. 
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4.8-9.1 GHz, 22.4 GHz, and 36.5 GHz, respectively using the 100-m Effelsberg telescope. 

5 Bell et al. (2000) estimated r] B = 0.4 ± 0.04 and t]a = 0.3 ± 0.03 at a frequencies of 6.0 and 17.6 GHz, respectively using the 
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3. Reduction and Analysis 

Before any profile fitting was done, the continuum level was subtracted from each spectrum, using up 
to a third-order polynomial fit to the baseline. This removed any instrumental slopes in the bandpass so 
that the baselines would be flat enough for a profile analysis, yet the shape of the original RRL would be 
preserved. All RRLs were processed with the NRAO GBTIDL data reduction package according to standard 
reduction protocol (see, e.g., Garwood et al. 2007). Particular care was taken to preserve the shape of the 
line wings. 

The equations used to generate Gaussian and Lorentzian profile fits for the RRLs in this survey were 

4>g(v) = A exp |- (V - j, (1) 

where the Gaussian FWHM is given as 

FWHMg(ct) = 2a V21n(2) (2) 

and 



' (y - v ) 2 + y 2 



Mv)=A - / (3) 



where the Lorentzian FWHM is given as 



FWHM L (y) = 2y, 



(4) 



where Ao represents line intensity, vo represents frequency at line center, v is the varying frequency across 
the breadth of the RRL, cr is the one-half value at lie intensity of the Gaussian line profile, and y is one-half 
value at 1/2 intensity of the Lorentzian line profile. 

In general, astronomical environments produce RRLs that experience both Doppler and pressure broad- 
ening. However, depending on the RRL transition and the physical properties of the gas, one type of broad- 
ening mechanism is usually observed to dominate over the other. At radio frequencies, Doppler broadening 
is almost always measured. In cases for which it is possible to detect both Doppler and pressure broaden- 
ing, the observed spectral line profile is expected to be a convolution between the Gaussian and Lorentzian 
profiles. The convolved shape is referred to as a Voigt profile and may be described by the formula 
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where Erfc is the complementary error function. 

Each H RRL presented in this article was fitted with Gaussian, Lorentzian, and Voigt profiles produced 
by the Mathematica NonLinearRegress least-squares-fitting routine. The various fit parameters and 
error estimates for each RRL are summarized in Table 1. Gaussian, Lorentzian, and Voigt models were 
generated using Equations (1), (3), and (5), respectively. For the Gaussian and Lorentzian fits we list the 
peak intensity, the FWHM line width, and the sum-of-squares (SOS) residual value, which measures the 
deviation of each point from the best fit model. For the Voigt fit we list the peak intensity, the Gaussian and 
Lorentzian FWHM line widths, and the SOS residual value. Also listed is the ratio a = FWMH L /FWMH G , 
which is the ratio of the Lorentzian-to-Gaussian component FWHM values and illustrates the degree of 
deviation of the line profile from a pure Gaussian or pure Lorentzian line shape. The estimates of the model 
parameters were found from a least-squares fit, minimizing the SOS residuals, and the model parameters 
are reported over a 99.7% confidence interval, representative of a 3o~ error on the fit. These errors, reported 
in Table 1, are on the last significant figure. In every case, the Voigt fit gives the smallest SOS residual, 
suggesting that it may be the best fit profile for the observed data. 

Figure 1 contrasts an apparently Gaussian profile (the H69a transition at 19590 MHz) with a strongly 
Voigt-shaped profile (the H109o' transition at 5009 MHz). For the H69a line, the best fit Voigt profile is 
almost exactly coincident with the Gaussian profile, indicating that pressure broadening is not significant 
for this case. In contrast, for the H109ck profile, the line intensity and shape clearly require contributions 
from both pressure and Doppler broadening to match the data. Figures 2-4 show the RRL baseline wings 
with the various profile fits in order of decreasing n. The evolution in line shape is clear, showing the 
shift from predominantly Voigt-shaped profiles at lower frequencies (4744, 5009, 5444, and 8584 MHz) to 
more Gaussian-shaped profiles at higher frequencies. Physically, this corresponds to a transition from lines 
broadened by a combination of both Doppler motions and electron impacts to lines broadened primarily 
by Doppler motions alone. The wings of the highest frequency RRL transition (24510 MHz) appear from 
visual inspection to be well-modeled by a Gaussian profile. However, between the 8873 MHz and 19590 
MHz transitions, the RRL profiles appear to fluctuate between either having shapes mostly described by 
Gaussian profiles (14129, 15281 MHz) and having wing shapes and breadths intermediate to those expected 
for Gaussian and Voigt profiles (9173, 10161, 17992, and 19590 MHz). As a general trend, we see support 
of the theory that for higher values of n, pressure broadening becomes increasingly important. 

Further statistical analysis was carried out using custom code within the SAS environment to compare 
the observed data to the Voigt, Gaussian, and Lorentzian fits for each of the Una RRLs listed in Table 
1. Correlation analyses between the observed data and that from the three distributions at all frequencies 
indicate in every case a strong correlation between the observed data and the data generated from the fitted 
statistical distribution (p < 0.0001). The median test (Table 2) examines whether there is a significant 
difference between the observed data and the generated data by counting the number of values in the data 
set above and below the calculated median and statistically comparing the frequencies. If the p-value is less 
than 0.05, this implies that there is a statistically significant difference between the observed data and the 
data from the proposed distribution. If the p-value is greater than 0.05, on the other hand, this implies that 
the two data sets come from the same distribution. For all H RRL frequencies, the Voigt distribution shows 
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that there is no significant difference between the observed data and the Voigt fit. This is particularly the 
case with the lowest frequencies (4744, 5009, 5444, 8584, and 8873 MHz), where the test suggests that there 
is no evidence at all that the observed data and the Voigt data are not from the same distribution. In fact, 
the Voigt data are considered the best representation for the observed data at all frequencies except 10161 
MHz, where the Gaussian data are considered better. We note that the 10161 MHz line appears to suffer 
from considerable spectral contamination on its positive velocity wing and that in all cases, the presence of 
the corresponding Hea RRL interferes with the fit on the negative velocity side of the Ha RRL. The analysis 
also indicates that the Lorentzian distribution is a significantly poor fit (p < 0.05) for all frequencies except 
4744 MHz, where it is not as good as the Voigt distribution (p=0.3323 compared to p=0.99999) but that it 
is statistically better than the Gaussian distribution (p=0.0878). 

The median test, however, is less robust than the Wilcoxon Mann Whitney U test (Table 2), which 
not only considers the values that differ from the median but their also their magnitudes. Examination of 
the Wilcoxon statistics indicates that again the Voigt is the better choice for the lower frequencies (4744, 
5009, 5444, 8584 MHz) but not necessarily so for the higher frequencies. For the profiles at 8873 MHz and 
above, the Voigt data are significantly different from the observed data. The Gaussian data are considered 
statistically closer to the observed data for the 14129, 15281, and 24510 MHz profiles. However, it is 
apparent that the remaining frequencies (9173, 10161, 17922, and 19590 MHz) are not suitably represented 
by any of the three distributions as the p-value is always less than 0.05, suggesting that the observed data are 
significantly different from the distributions data for these frequencies. The analyses also confirm previous 
results that the Lorentzian distribution is a poor fit for all frequencies except 4744 MHz, where it is not as 
good as the Voigt (p=0.1211 compared to p=0.4924) but is statistically better than the Gaussian (p=0.0006). 



4. Discussion 

RRL theory predicts that emission from ionized gas will produce Voigt line profiles, which consist of 
a Doppler (Gaussian) component and a pressure broadened (Lorentzian) component (Griem 1967). This 
assumes that no large-scale motions such as expansion or rotation exist within the observed region. Pressure 
broadening is caused by impacts from electrons and depends on the local electron density and the principle 
quantum number n. Pressure broadening is a sensitive function of n; as n increases, the line profile is 
expected to go from being dominated primarily by the Gaussian component to being affected by a significant 
Lorentzian component as well (Brocklehurst & Seaton 1972). This trend can be detected in several ways. 
Since the Lorentzian profile decays as a power law, rather than exponentially as the Gaussian does, RRL 
profiles in regions of ionized gas for which pressure from electron collisions is significant are expected to 
have increasingly broadened wings toward their baselines (see Figures 2-4). As theorized by Griem (1967), 
the amount of pressure broadening present can be quantitatively measured by the ratio of Lorentzian-to- 
Gaussian influences on the spectral line's FWHM value (i.e. parameter a in Table 1). By inspecting these 
values, we see that it does indeed generally increase with n. That is to say, the Lorentzian component 
becomes significantly dominant, which allows the wings of the Voigt profile to be detected. That these 
particular RRLs are indeed Voigt in nature is confirmed through both median and Wilcoxon statistical tests. 
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Since the frequency of an RRL is a function of n, the telescope beam size (and possibly the opacity 
as well) change with n. Detailed models of the electron temperature and density are required for a proper 
physical analysis of RRLs such as those reported here. Such a rigorous test of pressure broadening on 
a source as complex as Sagittarius B2(N) is beyond the scope of this paper and might prove challenging 
with the current data set. Nevertheless, the expected theoretical trends of pressure broadening appear to be 
correct, and we conclude from visual and quantitative analyses that at least four of the RRLs detected in 
this survey-4744, 5009, 5444, and 8584 MHz-are verifiably Voigt in nature. With today's spectrometers, 
sacrificing velocity resolution for more channels (leading to poor velocity resolution) is no longer an issue 
against detecting extended wings in line profiles. Spectrometer improvements to existing (e.g., EVLA) 
and future (e.g., ASKAP) interferometers hold much promise as probes of pressure broadening within high 
density, compact HII regions, and data achievable with the GBT will certainly provide a crucial role towards 
advancing our understanding of the larger-scale, extended ionized structures in the ISM. 
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Table 1 : Parameters of the best fit Gaussian, Lorentzian, and Voigt profiles of the selected RRLs. Note - v represents the approximate RRL frequency in MHz, 
and n represents the principal quantum number of the RRL observed. A v , Ac, and A L represent the line intensities of the Voigt, Gaussian, and Lorentzian profiles, 
FWHM(G) and FWHM(L) represent the respective contributions of the Gaussian and Lorentzian line widths to the Voigt line width, a is the ratio of line widths: 
FWHM(L)/FWHM(G). SOS(V), SOS(G), and SOS(L) are the sum-of-squares values for the Voigt, Gaussian, and Lorentzian profiles. FWHM C and FWHM L represent 
line intensity and FWHM values for the Gaussian and Lorentzian profiles. Values given in parentheses are the 3cr confidence intervals of the fits. 
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Table 2: Fit quality statistics for various RRL model types. Note-Strong correlation (p < 0.0001) was 
measured between each observed RRL and its corresponding statistical fit. According to the median test, the 
Voigt function produces the best fit profile in every instance except for the RRL at 10161 MHz. According 
to the Wilcoxon test, the RRLS at 4744, 5009, 5444, and 8584 MHz are best described by the Voigt profile, 
the RRLs at 14129, 15281 and 24510 MHz are best described by the Gaussian profile, and the RRLs at 
9173, 10161, 17922, and 19590 MHz are not well-described by any of the three distributions investigated. 
See S3 for details. 




Fig. 1. — The observed H 109a RRL at 5009 MHz (right) and H 69a- RRL at 19590 MHz (left) from observations taken towards Sagittarius B2(N) illustrating the 
shift from nearly a pure Doppler (Gaussian) line profile to one with a substantial amount of pressure broadening (Lorentzian). In both cases, the Voigt profile has a lower 
SOS residual value than either the Gaussian or the Lorentzian profile. The observed data are shown in the bold trace. The best Voigt fit is shown as a solid line, the best 
Gaussian fit is shown as a dotted line, and the best Lorentzian fit is shown as a dashed line. The weaker line to the negative velocity side of the H RRL is the corresponding 
He RRL. 
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Fig. 2. — Wings of RRLs between 4744 MHz and 8584 MHz. View of full raw data is inset (y-axis ranges from to 1.5 K, x-axis 
ranges from -200 to 200 km s~'). The observed data are shown in the bold trace. The best Voigt fit is shown as a solid line, the best 
Gaussian fit is shown as a dotted line, and the best Lorentzian fit is shown as a dashed line. In all four cases, the best statistical fit 
to the data is the Voigt profile. 
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Fig. 3. — Wings of RRLs between 8873 MHz and 14129 MHz. View of full raw data is inset (y-axis ranges from to 1.5 K, 
x-axis ranges from -200 to 200 km s" 1 ). Figure labels, see Figure 2. 
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Fig. 4. — Wings of RRLs between 15281 MHz and 24510 MHz. View of full raw data is inset (y-axis ranges from to 2.0 K, 
x-axis ranges from -150 to 150 km s" 1 ). Figure labels, see Figure 2. 



-15- 



REFERENCES 

Balser, D. S., Bania, T. M., Rood, R. T. & Wilson, T. L. 1999, ApJ, 510, 759 

Bania, T. M., Balser, D. S., Rood, R. T., Wilson, T. L. & Wilson, T. J. 1997, ApJS, 113, 353 

Bell, M. G., Avery, L. W., Seaquist, E. R. & Vallee, J.P. 2000, PASP, 112, 1236 

Brocklehurst, M. & Seaton, M. J. 1972, MNRAS, 157, 179 

Brown, R. L., Lockman, F. J. & Knapp, G. R. 1978, ARA&A, 16, 445 

Chaisson, E. J. 1973, ApJ, 186, 555 

Churchwell, E. 1971, A&A, 15, 90 

Churchwell, E. & Walmsley, C. M. 1975, A&A, 38, 451 

Foster, T. J., Kothes, R., Kerton, C. R. & Arvidson, K. 2007, ApJ, 667, 248 

Garwood, R. W., Braatz, J. A., Maddalena, R. J., Marganian, P., & Minter, A. H. 2007, ASPC, 376, 277 

Gordon, M. A., Sorochenko, R. L. 2002, Radio Recombination Lines: Their Physics and Astronomical 
Applications (Kluwer) 

Griem, H. R. 1967, ApJ, 148, 547 

Kantharia, N. G., Anantharamaiah, K. R. & Payne, H. E. 1998, ApJ, 506, 758 
Keto, E., Zhang, Q. & Kurtz, S. 2008, ApJ, 672, 423 

Mehringer, D. M., Palmer, P., Goss, W. M. & Yusef-Zadeh, F. 1993, ApJ, 412, 684 
Miao, Y., Mehringer, D. M., Kuan, Y.-J. & Snyder, L. E. 1995, ApJ, 445, L59 

Qin, S.-L., Zhao, J.-H, Moran, J. M., Marrone, D. P., Patel, N. A., Wang, J. -J., Liu, S.-Y & Kuan, Y-J. 2008, 
ApJ, 677, 353 

Quireza, C. Rood, R. T., Bania, T. M., Balser, D. S. & Maciel, W J. 2006, ApJ, 653, 1226 

Reid, M. J., Schneps, M. H., Moran, J. M. & Gwinn, C. R., 1988, ApJ, 330, 809 

Remijan, A. J., et al. 2007, BAAS, 39, 963 

Roelfsema, P. R., & Goss, W. M. 1992, A&A Rev., 4, 16 

Rohlfs, K. & Wilson, T. L. 1996, Tools of Radio Astronomy (Springer- Verlag) 

Shaver, P. A. 1980, A&A, 91, 279 



-16- 



Shaver, P. A., McGee, R. X., Newton, L. M ., Danks, A. C. & Pottasch, S. R. 1983, MNRAS, 204, 53 
Simpson, J. P. 1973, Ap&SS, 20, 187 

Smirnov, G. T., Sorochenko, R. L. & Pankonin, V. 1984, A&A, 135, 116 
Snyder, L., Kuan, Y.-J., Miao, Y., Lovas, F. 1995, ASPC, 74, 107 

Stepkin, S. V., Konovalenko, A. A., Katharia, N. G. & Shankar, N. U. 2007, MNRAS, 374, 852 

Ulich, B. L., Haas, R. W. 1976, ApJ, 30, 247 

Wink, J. E., Wilson, T. L. & Bieging, J.H. 1983, A&A, 127, 211 



This preprint was prepared with the AAS MgX macros v5.2. 



